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Abstract 



A simple solution to the BFKL equation is obtained as a series in the number 
of real gluons emitted with transverse momentum greater than some small cutoff 
/i. This solution reveals physics inside the BFKL ladder which is hidden in the 
standard inclusive solution, and lends itself to a straightforward Monte Carlo im- 
plementation. With this approach one can explore new useful physical observables, 
which are shown to be independent of the cutoff fi. In addition, this approach al- 
lows the imposition of kinematic constraints (such as energy conservation) which 
are important at finite energies. The distribution of J2 E± of particles in a central 
rapidity bin between two widely-spaced jets is presented as an example. 



The BFKL (Balitsky-Fadin-Kuraev-Lipatov) equation [|T| , which systematically 
resums powers of as times large rapidity intervals or logarithms of Feynman x in 
perturbative QCD, has recently moved from the purely theoretical realm to the 
phenomenological arena. The classic prediction of the BFKL resummation is the 
rise of ^2(0;) at small x in deep inelastic scattering. Unfortunately, due to the 
resiliency of the parton density functions and DGLAP evolution equations, this 
observation |^, ^ is still open to interpretation Alternatively, by tagging jets 
on both ends of a large rapidity interval, as suggested by Mueller and Navalet 
0, it is possible to unambiguously isolate the effects of the BFKL ladder from 
the parton density functions. Predictions of this kind include the rise in the cross 
section as a function of the rapidity interval, both in hadron-hadron |^ and lepton- 
hadron colliders 0. In particular, recent preliminary results from HI |^ on deep 
inelastic scattering with a tagged jet at fixed Xj show an intriguing rise in the cross 
section with decreasing xbj that is not well-explained either by a fixed-order matrix 
element calculation |^ or by parton shower Monte Carlo simulations. The BFKL 
prediction however, seems to be in good agreement. Other observables related 
to the kinematics of the tagging jets, such as the decorrelation in azimuthal angle 
Tod , have also been considered. 



With the advent of BFKL phenomenology, it has become a necessity to under- 
stand the range of validity and the errors inherent in the BFKL approximations. 
The calculation of the next-to-leading logarithmic (NLL) corrections to the BFKL 



matrix elements is currently under way [[T^ . Short of these full NLL corrections, it 



has been shown in [jT2[ that some of the largest corrections to the asymptotic the- 
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ory are purely kinematic in origin. In the standard solution to the BFKL equation 
the transverse momenta of the ladder gluons are integrated from zero to infinity. 
However, at physical energies the inclusion of kinematic constraints on these inte- 
grals can significantly modify the predictions of the theory, even if the difference is 
formally subleading in the asymptotic expansion. In Ref. |jl2| an effective rapidity 



interval was defined as an attempt to include some of these kinematic effects in the 
BFKL calculation. 

In this paper we present a new solution to the BFKL equation in a form natu- 
rally suited for physical interpretation and Monte Carlo implementation . This 
approach also offers a natural way to impose kinematic constraints and to assess 
the uncertainties due to the asymptotic nature of the equation. The key ingredient 
to this solution is the introduction of a lower cutoff fi on the transverse momentum 
of the real gluons that are produced. However, the cutoff is introduced in such a 
way that, for sufficiently small /i, any infrared safe observable is independent of the 
cutoff. In the inclusive case the Monte Carlo method generates an exact answer, 
identical to the known solution in the literature, with no arbitrary parameters. In 
addition, new experimental observables which depend on the observation of the 
ladder gluons can be obtained |T^ . 



Let us begin with a general description of a semi-hard process, where the parton- 
parton center-of-momentum energy is much larger than the typical momentum 
transfer scale Q. In this limit the partonic cross section factorizes into the form 

da 



K(Pa±) fiyab,Pa±,Pb±) 
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A physical interpretation of this form of the cross section is represented in Fig. 1. 
The process consists of two distinct scatterings, which occur at widely-separated 
rapidities, ?/„ and i/b, and small transverse momenta Pa± = \Pa±\ ~ Pb± = IPbJ-l ~ Q- 
Each scattering form factor, K(Pa±) and Vb{pb±), depends only on the transverse 
momentum which flows into its particular vertex. The precise form of the form 
factors, however, depends on the specific partons involved in the scatterings. The 
function f {Uab, Pa±, Pb±) , which connects the two scatterings, is essentially a prop- 
agator which allows pb± to flow to Pa± by emitting gluons over a rapidity interval 
Vab = Ua — Ub ^ ^^{s / Pa±Pb±) ■ This fuuctiou is uuivcrsal, and it is our object of 
interest . 

The form of the cross section given in (|I]) is valid in the limit of large t/ab- 
In that limit the perturbative contributions to f {Vab, Pa±, Pb±) th.a.t are leading in 
asVab can be resummed systematically with the aid of the BFKL equation. We now 
present the equation and describe the physics that it encodes. It can be written 

where as = c^sNc/tt- The boundary condition for the equation is 

f{0,Pa±,Pb±) = ls^'\Pa± + Pb±) , (3) 

which corresponds to no gluon emissions and enforces conservation of transverse 
momentum. The flrst term on the right-hand side of the BFKL equation (^, upon 
iterating from the boundary condition n times, gives the squared amplitude for 



producing n real gluons, in the approximation that the gluons are well-separated in 
rapidity. The amplitude with three gluon emissions is represented by the Feynman 
diagram in Fig. 2. The second term on the right hand side of (|^) gives the virtual 
corrections which "reggeize" the t-channel gluon propagators, represented by the 
heavy solid lines in the figure. 

An important point here is that the singularities in the real and virtual terms 
of the integrand cancel as k± becomes small. Thus, we can cut off the integral at 
k± = fi for both the real and virtual gluons, and the neglected contribution will 
vanish as /i — > 0. With this cutoff we can explicitly do the integration over the 
virtual gluon corrections, obtaining: 



where the integral over real gluons is restricted to k± > fi. From hereon we neglect 
the terms of 0{fi^ /p^jj. Then we can simplify this equation further by making the 
substitution 



with the same boundary condition (^. 

In we see that the dependence of the virtual corrections on fi simply ex- 
ponentiates into an overall factor. Starting with the boundary condition, we can 
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iterate equation (H) to obtain a series solution / = I]^o where the nth term is 
the contribution from the emission of n real gluons having k±_ > fi. This contribu- 
tion can be written as a product of integrals over the phase space of each gluon in 
the ladder: 



riyab,Pa±,Pb±) 



n 

1=1 



otsdy, 



/i 



2 \ O'sVi+la 



(7) 



X -5^'^\pa± + qr^±) 



with 



qj± = Pb± + J2^i 



i=l 



and qo± = Pb±- The integrals in rapidity are ordered with yi, = i/q < yi < 
■ ■ ■ < yn < Vn+i = Ua- Written in this manner, the solution recovers the simple 
interpretation as the emission of real gluons by the exchange of "reggeized" gluons 
in the t-channel as in Fig. 2, with /i the infrared cutoff to the gluon Regge trajectory. 
For clarity, we emphasize that only the gluons exchanged in the t-channel are 
"reggeized" , while the emitted partons are standard gluons. 

Each term in the series is positive definite. Therefore, it is straightforward to 
implement this solution as a Monte Carlo simulation. Given pbj_ and the rapidity 
interval ya — Vb, we first sample the distribution in the number n of gluons in the 
ladder. Next, we produce the four-momenta of the n gluons successively as given by 
the distribution (|^). Finally, we fix pa± by conservation of transverse momentum. 
In practice the events are produced using approximate distributions and are then 
reweighted. 
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To understand more clearly the cutoff dependence of this solution, let us cal- 
culate the quantity F = J f dp'^j_d(f)a, where the integral just fixes Pa± via the 
(5-function in (|^). For simplicity of the discussion here, it is convenient to consider 
a modified equation, obtained by replacing qi± pb± everywhere in (0) and by 
setting the upper limit on the integrations to k'j_ = pl^. Now the nested integrals 
can be done analytically, the series can be summed, and the dependence on fi van- 
ishes identically. In fact we obtain F = 1, and the distribution in the number of 
ladder gluons is just a Poisson distribution with mean 

(n) = asyab^T^ipl±/lJ''^) ■ (9) 

The physical significance of the cutoff is now apparent. As we lower /i, the number 
of gluons emitted in the rapidity interval grows logarithmically. However, the 
average k± of each gluon is reduced in such a way that, for suitable infrared-finite 
observables, the dependence on the cutoff vanishes. In the exact solution (0) the 
distribution in the number of ladder gluons will differ somewhat, but the qualitative 
features of this discussion still apply. 

We now present in Fig. 3 a plot of the Monte Carlo solution compared to the 
standard BFKL solution 

1 °° poo 

f{yab,Pa±,Pb±) = TT^-^ E e^"(^---)/ rfz/e^("'^)^-(pL/pL)^% 

(10) 

with (t)ab = <Pa-<Pb and 

Inl + 1 

a;(n,z/) = 2a,[^(l)-Re^(^— + ZZ/)] , (11) 
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where ip is the logarithmic derivative of the Gamma function. In this figure we fix 
Pb± = 50 GeV and i/ab = 4, and we plot J f dcp as a function of pa±. The Monte 
Carlo solution agrees with the standard solution ([10|) for this curve as long as /i is 
smaller than the bin size used near the peak. 

Next, we calculate an observable for the Tevatron at ^/s = 1800 GeV that 
exhibits the full potential of the Monte Carlo approach to BFKL physics. We tag 
on the two jets a and b with the largest and smallest rapidities with Pa±,Pb± > 20 
GeV. Then we sum the E±_ of all the particles within a bin in rapidity and plot 
the distribution in = J2E±. In Fig. 4 we plot the distribution of summed 
in the bin \y\ < 0.5 for fixed values of i/a = 2.5 and yi, = —2.5 for the tagging jets. 
In this calculation we have also made an improvement to the BFKL prediction 
by including the kinematic contribution of all of the physically-produced particles 
to the Feynman-x parameters in the parton density functions of the proton and 
anti-proton. In practice, this makes a large effect due to the constraint of total 
energy conservation imposed by the parton density functions. More details of this 
kinematic improvement and other phenomenological results will be reported in an 
expanded paper . 



In conclusion, we have presented a solution to the BFKL equation for large 
rapidity intervals as a series in the number of real gluons emitted above a cutoff 
/i. We have incorporated this solution into a Monte Carlo simulation, and we 
have shown that it reproduces exactly the BFKL dynamics with no dependence on 
the cutoff for small fi. An advantage of a Monte Carlo solution to BFKL is that it 
allows one to study the effects of the gluons emitted in the middle of the ladder and 
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even to make experimental cuts on these ladder gluons. Finally, with the Monte 
Carlo simulation it is also possible to improve the convergence of the resummation 
by including kinematical effects exactly in the cross section. 

The author would like to thank Mike Albrow, Vittorio Del Duca, Terry Heuring, 
Joey Huston, Harry Weerts, and C.-P. Yuan for useful conversations and Wu-Ki 
Tung for a critical reading of this manuscript. 
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Fig. 1: A schematic picture of the cross section for producing particles at large 
rapidity separation. 



Vb < yi < y2 < VS < Va 




Fig. 2: The amplitude for three gluon emissions in the BFKL approximation. 
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Fig. 3: Comparison of the Monte Carlo solution (histogram) with the standard 
solution to BFKL (solid). 
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Fig. 4: The distribution of Sj^ as explained in the text. 
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